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Abstract 

In the turbulent boundary layer above a flat plate, the velocity profile is known 
to have the form v = Vo[-lnz + constant]. The distance from the wall in dimen- 
sionless units is z and vq is an uniquely defined velocity scale. The number k is 
universal and measurements over several decades have shown that it is nearly 0.42. 
We use a randomly stirred model of turbulence to derive the above law and find 
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~ 0.52. 
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The turbulent flow over a stationary plate has been investigated in great detail for 
well over a century. One of the important quantities studied is the velocity profile and in 
particular a good deal of attention has been paid to the velocity in the boundary layer. 
The existence of the boundary layer was first postulated by Prandtl [1,2]. Due to the 
presence of viscosity, the flow speed has to be zero on the surface of the plate. Within 
a short distance, the velocity has to rise to its bulk value. This small distance in which 
the velocity experiences a rapid change is called the boundary layer. Prandtl observed 
that in the boundary layer the velocity v scales with the distance z from the boundary. 
Adjoining the boundary is an extremely thin laminar boundary layer, after which comes 
a buffer zone and then the turbulent boundary layer. In this regime one has [2,3] 

- = -[ln{z/z Q )i (1) 

a result known as the law of the wall, where v o is a reference velocity to be defined below 
and k is an universal constant which so far has only been experimentally measured [4]. 

With the coming of Kolmogorov's theory [5] of homogeneous isotropic turbulence, 
where simple dimensional argument yielded remarkable results, dimensional arguments 
were similarly used in this problem to obtain eqn.(l). A more careful analysis [6,7] of 
the Navier Stokes equation also supported the above result, which is one of the corner- 
stones in the theory of the turbulent boundary layer. Over the last two decades, a lot 
of effort has been spent in providing a more microscopic justification for Kolmogorov's 
scaling in homogeneous isotropic turbulence, thereby enhancing our understanding of the 
energy cascade involved. Analysis of homogeneous isotropic turbulence has thus become 
a problem of non-equilibrium statistical physics [8-14]. In this work, we propose to use 
techniques from statistical physics to study the law of the wall and in the process provide 
calculations for the universal constant k. 

We begin by explaining the geometry and the traditional formulation. The flat plate 
(infinitely large) sits at z = (in the x-y plane) and the flow (in the steady laminar 
situation) occurs in the x-direction, due to some maintained pressure gradient in that 
direction. We can view this as a two dimensional problem - the dynamics occurs only in 
the x-z plane. Continuity (for an incompressible flow) implies + ^ = 0. Translational 
invariance in the x-direction makes quantities x-independent and hence |j = 0. The 
resulting constant value of w is found to be zero since w = at z = 0. In the laminar 
regime (low Reynold's number) the velocity u satisfies y 2 u = and since the only 
dependence of w is on z, we have = or a linear profile for u, i.e. u = uo~, where 
Uq and a are constants. This is the usual Couette flow. 

We now pass to high Reynold's number where the non-linear terms in the Navier 
Stokes' equation need to be retained. We assume that the flow is turbulent, so that it 
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pays to write the velocity field as v = V + u, where < v > = V. This has the implication 
that u (components u x ,u y ,u z ; which from now onwards we will designate as Ui,u 2 ,U3 
respectively) is a purely fluctuating component and V (components V x and V z , henceforth 
to be called as V\, V3 respectively) is the mean flow. No mean quantity will depend on x 
due to translational invariance and hence V z — i.e. V3 = 0, by the same argument as 
given above. Hence finding the flow profile is the same as finding V\{x^) (i.e. V x (z)). 
Writing the Navier Stokes' equation as 

9v VP 2 - 

+ ( v .V)v = + v\7 2 v, (2) 

01 p 

substituting v — V + u and taking an average yields (for very high Reynold's number) 
after standard manipulations [7] 

v-— = < U1U3 > + constant (3) 
oz 

At this point, one invokes dimensional arguments to argue V\ = ^ln(z/ ' zq). Our goal is 
to calculate < u±u 3 > directly. 

We proceed by writing an equation of motion for the fluctuating velocity field u. This 
is the usual Navier Stokes' equation supplemented by a random forcing term which is 
generated by the interaction of the mean flow and the turbulent flow and the boundary 
conditions. We write this as 

dui d 1 dP 2 „ , AS 

^ +U ^ = - p d X - + US7 + (4) 

with = (this determines P in terms of u { and makes a non-linear term), where 
fi is the random force, whose correlation is to be specified. For the Kolmogorov theory 
of homogeneous isotropic turbulence, this forcing gives rise to spatially correlated noise 
[9-14] and in a D-dimensional space, the Fourier components fi(k,u) of fi(x,t) satisfy 
< fi(k 1 ,oj 1 )f j (k 2 ,u;2) >= F(k l )5(k 1 +k 2 )5(uj l +Lo 2 )P ij (k l ) where F(k) ~ £r D and i^j(fc) is 
the projection operator. In co-ordinate space < fi(xi, ti)fj(x 2 , t 2 ) >oc PijF{x\ 2 )5{ti—t 2 ) 
where F{x\ 2 ) is a constant for x\ 2 — > 0. The combination of the non-linear terms and 
the forcing changes the basic relaxation rate in the problem. In momentum space, in 
the absence of the non-linear terms, the basic frequency is clearly uk 2 . The non-linear 
terms constitute a relevant perturbation and makes v diverge for small /c.For the forcing 
of k~ D , one gets the Kolmogorov time scale where v = T^k~ 4 ^ 3 and the basic frequency is 
TTq/c 2 / 3 when the sweeping effects are removed. For this problem of calculating < uiu 3 >, 
sweeping effects will set the time scale (interaction of large and small eddies) and hence 
the basic frequency scale is u Q k [13,14]. It should be noted that for homogeneous isotropic 
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turbulence the cross correlation < ui(f,t)u 3 (f,t) > will be zero - the finite value of the 
correlation is because of the breaking of isotropy by the rigid boundary at z — 0. 

To deal with the situation, restricting ourselves to z > 0, we need to use field variables 
for which we introduce Fourier transforms in the x-y plane and keep the co-ordinate 
description along the z-axis. To take into account the solenoidal condition, it is convenient 
to work in the potential A, so that u — V x A. We expand A as 

A(r,t) = ]T Mk,z,uj)e^-^ (5) 

k,u> 

and thereafter write down the equation of motion for A from eqn.(4). We will denote the 

— * — * — * 

random force in the equation of motion for A by F. Specifying the correlations of F will 
now be the main task. One part of the correlations of Fi and Fj will be the homogeneous 
isotropic part, giving < Fj(&, zi)Fj(—k, z 2 ) > = 5ij-^F(z 12 ) where F(z 12 ) is a constant for 
the spatially correlated noise. This reproduces the Kolmogorov spectrum. The boundary 
at z — 0, will cause anisotropy in the correlation of Fj in the x-y plane (this anisotropy 
causes the correlation < U\U 3 > to be non-zero in the presence of the wall). Considering 

— * 

all these factors, the final correlations of F assume the form 

< F i {k,z l )F j {k,z 2 ) > = A-j^ [8ij + 5n5 j3 + 5 i2 5 j3 \ (6) 

where e is the rate of production of turbulent energy. The first term on the right hand 
side produces the isotropic Kolmogorov response for the fluctuating field. The second and 
third terms reproduce the effect of the boundary. 

We now start with a linear equation of motion for Ai(k,z) with the effect of the 
non-linear terms included in an effective viscosity v e ff- This reads 

+ „„,(*» = * (7) 

with the correlations of Fj specified in eqn.(6). We see from eqn.(3) that our primary goal 
is to calculate the stress tensor t\z = < u\u 3 >. This can be written as 

f d 2 k duj , 
< ui(r,f)u 3 (r,f) > = J j^^— <u 1 (k,z,u))u 3 (-k ) z ) -u))> 



d 2 k duo , 2 . . . . 

k 2 < A 3 {k, z, u)Ai(— k, z, — uo) > 



= -Ae j M [ G ^l G -^ d z 1 dz 2 

J (2tt) 2 2tt A; 4 J v 2 ef} 

(8) 

where G±(z, zi) are the Green's functions associated with the operator L = — k 2 =F 
where fi = uj/v. In arriving at the last line of eqn.(8), we have used the correlation of 
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eqn.(6) and noted the fact that the angular average of k\ is zero. Defining k± = k 2 ± ifl 
and noting that we are interested in a single rigid boundary, G± can be written down as 
[17] 

G ± (z,z,) = smh ^<) ^ ±Z> (9 ) 

K± 

Carrying out the integration in eqn.(8), we have 

d 2 k du k\ 1 



< M1M3 > = — Ae j 



(2tt) 2 2tt A; 4 k\k 2 _ 
2 — (1 - e- K + z )(l - e- K - z ) (10) 



U eff 



For large z, the stress tensor tends to a constant value — V 2 which is given by 

d 2 k du cos 2 9 1 1 



r d*k 

1 W) 

I 



V ° Ae I (2tt) 2 2tt k 2 u 2 eff k\k 2 

At r d 2 k cos 2 e 1 



(2n) 2 u eff k* 



[The integrals have to be cut-off at a low momentum ko] 

In the sweeping time dominated situation u e ff(k) = VQe 1 ^k" 1 k Q and the above integral 
needs to be cut-off at a lower limit ko leading to Vq = It is this Vq which sets 







the scale for the velocity in the boundary layer. We now need to specify the region of 
z-space in which we are interested. As already discussed, the law of the wall in eqn.(l) 
holds somewhat away from the boundary. The relevant length scale is k^ 1 = = 

N/87r i/ o //2 ^i72 ■ O ur interest is in the value of z for which kr,z > 1 and from eqn.(10), we 
find the correction to V 2 is 

(12) 

where z + = koz (Note that for small k we can set e~ k ° z ~ 1). The quantity A/u 3 is 
a universal number [10-12] known from the study of isotropic homogeneous turbulence. 
Using that universal value we find 

< Ul u 3 >= -V 2 [l-l^-^] (13) 

We now return to eqn.(3) and observe that the constant can be calculated by going far 
away from the wall, where the viscous stress is negligible. This leads to 

v— = < Ul u 3 > + V 2 (14) 
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(A more careful calculation would yield a correction term in the above equation which 
vanishes for large Reynold's number). Using eqn.(13) for the fluctuating stress tensor and 
integrating, we get 



The experimental values of k cluster around 0.42. 

In summary, we have used a randomly forced Navier Stokes' equation to calculate the 
fluctuating stress tensor and have thereby obtained the logarithmic law of the wall. The 
strong point of the procedure is the ability to calculate the universal amplitude appearing 
in the velocity profile. Our calculated number is in reasonably good accord with the 
experimentally measured value. 
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